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Abstract. The fractional Laplacian can be obtained as a Dirichlct-to-Ncumann map via an ex- 
tension problem to the upper half space. In this paper we prove the same type of characterization 
for the fractional powers of second order partial differential operators in some class. We also got 
a Poisson formula and a system of Cauchy-Riemann equations for the extension. The method is 
applied to the fractional harmonic oscillator H" = (—A + |a:|^)'' to deduce a Harnack's inequality. 
A pointwise formula for H" f{x) and some maximum and comparison principles are derived. 



1. Introduction 

In the last years there has been a growing interest in the study of nonlinear problems involving 
fractional powers of the Laplace operator (— A)°', < cr < 1. The fractional Laplacian of a function 
/ : M" R is defined via Fourier transform as 

(1-1) (^/(e) = ier^/(0, 

and it can be expressed by the pointwise formula 

(1-2) (-Ar/(a:)=c„,.P.V. / M^I^dz, 

where Cn^a is a positive constant. Observe from (jl.2p that the fractional Laplacian is a nonlocal 
operator. This fact does not allow to apply local PDE techniques to treat nonlinear problems for 
(— A)'^. To overcome this difficulty, L. Caffarelli and L. Silvestre showed in [5] that any fractional 
power of the Laplacian can be determined as an operator that maps a Dirichlet boundary condition 
to a Neumann-type condition via an extension problem. To be more precise, consider the function 
u = u{x,y) : R" x [0, oo) R that solves the boundary value problem 

(1.3) u{x,0)^f{x), zeR", 

(1.4) Ar,U + ^ ^^'^ Uy + Uyy = 0, X S R" , ?/ > . 

Then, up to a multiplicative constant depending only on cr, 

- \im y'~^^Uy{x,y) = {-Arf{x). 

This characterization of (— A)'^/ via the local (degenerate) PDE (jl.4p was used for the first time in 
[T] to get regularity estimates for the obstacle problem for the fractional Laplacian. 

To solve (|1.3p - (|1.4p . Caffarelli and Silvestre noted that (|1.4I) can be though as the harmonic exten- 
sion of / in 2 — 2cr dimensions more (see [2]). From there, they established the fundamental solution 
and, using a conjugate equation, a Poisson formula for u. Furthermore, taking advantage of the 
general theory of degenerate elliptic equations developed by Fabes, Jerison, Kenig and Serapioni in 
1982-83, they proved Harnack's estimates for u (and thus for /). 
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Let fl be an open subset of M", n > 1, and let dr] be a positive measure defined on fl. Consider a 
linear second order partial differential operator L, that we assume to be nonnegative, densely defined, 
and self-adjoint in L'^{^1, drf). The fractional powers < ct < 1, can be defined in a spectral way, 
see Section [2] 

The aim of this paper is to describe any fractional power L'^ as an operator that maps a Dirichlet 
condition to a Neumann- type condition via an extension problem as in [2], developing also the cor- 
responding properties (Poisson formula, fundamental solution, conjugate equation, Cauchy-Riemann 
equations). With this characterization, the interior Harnack's inequality for any fractional power of 
one of the most basic Schrodinger operators, the harmonic oscillator H = —A -I- , is consequently 
deduced. Besides, we find an explicit pointwise expression for the nonlocal operator H"^ that will 
allow us to get some maximum and comparison principles. 

Fractional operators appear in physics, when considering fractional kinetics and anomalous trans- 
port [H] . 

The heat-diffusion semigroup *^}oq generated by L will play a crucial role in our work. 
Our first main result is the following. 

Theorem 1.1. Let f G Dom(L'^). A solution of the extension problem 

(1.5) u{x,0)^f{x), onQ; 

(1.6) ^LxU H Uy + Uyy — 0, in Q X (0, oo); 

is given by 

(1-7) u{x,y) e-*^(LV)(a:)e-fe 



and 

, „x m(x, u) — u(x, 0) r(— cr) ^- , 1 1 9„ , N 

Moreover, the following Poisson formula for u holds: 

(1.9) u{x, y) ^ , , / e /(a;)e « -;— = / e /(a;)e — . 

All identities in Theorem 1 1.1 1 are understood in L^(0, dr]). Note that a solution u to the degenerate 
boundary value problem (|1.5p - (|1.6l) is written explicitly in terms of the heat semigroup e~*^ acting 
on L'^ f. From here, the Poisson formula (jl.9D can be immediately obtained (see the proof in Section 
[2]), where no fractional power of L is involved. When L = — A, the extension result of [2] is recovered 



(see Examples 12. 14p . More properties concerning the Poisson formula are contained in Theorem 12. II 
Moreover, (|1.9I) can be derived as in [5] (^Remark 12. 6^ : use the fundamental solution (that involves 
the kernel of the heat semigroup generated by L) and an appropriate conjugate equation (|2.10p to 
infer the Poisson kernel (see (|2.1ip ). The conjugate equation will be studied in detail by defining 
Cauchy-Riemann equations (I2.12p adapted to equation (|1.6p . See Section [5] 

If L has discrete spectrum, i.e. Lipk = Xk4>k, > 0, and {<pk}k£fiQ is an orthonormal basis of 
L^(f2, dri), the definition of the fractional power L'^ is given in the natural way: if / € L^{^, drf) has 
the property that Y.k ^l" I (/. <^fc) |^ = Efc ^l" |L /-^fc ^'7| < oo, then 



(1.10) L"f^Y.^k{f,'f>k)cl^k, sumini2(f7). 

k 



In Section [3] it is shown that, under this assumption, (jl.5p - (ll.6p has a unique solution u (vanishing 
as y oo) such that (|1.8p holds in the L^(ri)-sense. The proof is elementary using orthogonal 
expansions: just write u{x,y) = J2k'^k{y)(t>k{x) and observe that the coefficients Ck satisfy a Bessel 
equation. Hence, for the existence and uniqueness in this case, the general theory of degenerate PDE's 
mentioned above is not needed. This method also gives us local Neumann solutions (see Subsection 



EXTENSION PROBLEM AND HARNACK'S INEQUALITY 



3 



Let us now turn to the case of the fractional harmonic oscillator. We will be able to define H"^ f for 
all tempered distributions /. If / is a function that has also some local regularity then the extension 
result is true in the classical sense (Theorem 14.21 and Remark 14. 3p . This last fact is an essential 
ingredient for the second main result of this article: the interior Harnack's inequality for H"^ . 

Theorem 1.2. Let xq G M" and R > 0. Then there exists a positive constant C depending only on 
n, a, xq and R such that 

sup f < C inf /, 

for all nonnegative functions f : M" — > R that are in Bji{xq) and such that H'^ f{x) — for all 

X e Br{xo). 

The Harnack's inequality is valid for < a < 1 and the proof given in Section!?] is based (as we 
already remarked) on Theorem 14.21 and the Harnack's inequality for degenerate Schrodinger operators 
proved by C. E. Gutierrez in [3] (this idea is contained in 2, for the case of the fractional Laplacian). 
The Harnack's inequality ioi H (cr = 1) follows from general results (see [13]). 

The final part of the paper is devoted to the study of the pointwise expression of the fractional 
harmonic oscillator and some of its consequences. To that end we collect some previous facts about 
the fractional Laplacian (—A)"'. The natural way to arrive to (|1.2p starting from (|l.ip would be by 
taking the inverse Fourier transform. However, this path can be avoided if we consider the classical 
formula for L"^ that involves the heat-diffusion semigroup generated by L: 

1 r°° fif 

(1-11) ^'^(") = fF^X (^"'/(^)-/(^)) 

Note that (fLTT|) is motivated by the identity A'^ = /(,°°(e"*^ - 1) A > 0. When L = - A and 

/ e iS in (|l.lip . the Fourier transform recovers (|l.ip . Furthermore, the formula allows us to obtain 
(in a very simple way) expression (II. 2p with the constant Cn.a computed explicitly and in particular 
to see (Proposition 15. 3p that if a function / is around some x £ R" then 

lim (-A)^/(a;) = -A/(x). 

cr— Vl^ 

In Section [5] we put L = H in (jl.lip to get a pointwise formula for H°'f{x) (see Theorem 15. 7p and, 
from there, some maximum and comparison principles for . 

Throughout this paper S is the Schwartz class of rapidly decreasing C°°(R") functions, the letter C 
denotes a constant that may change in each occurrence and it will depend on the parameters involved 
(whenever it is necessary we point out this dependence with subscripts) and F stands for the Gamma 
function. We restrict our attention to < cr < 1 and, in this range, F(— cr) := < 0. 

2. The extension problem 

We begin with the basics of the spectral analysis that will be used throughout this Section. The 
complete details can be found in S] Ch. 12 and 13]. Since L is a nonnegative, densely defined and 
self-adjoint operator on L^[VL,drj) = L'^{n), there is a unique resolution E of the identity, supported 
on the spectrum of L (which is a subset of [0, oo)), such that 

A dE{X). 

The identity above is a shorthand notation that means 

/•oo 

{Lf,g)mn) = / A dEf,g{\), f e Dom(L), g G L\n), 
Jo 

where dEf^g{X) is a regular Borel complex measure of bounded variation concentrated on the spectrum 
of L, with d \Ef^g\ (0, oo) < ||/||i2(Q-) ||3|Il2(-q'). If 4'W is a real measurable function defined on [0, oo), 
then the operator (j){L) is given formally by 

(2.1) 0(i) = /"^(A) dE{X). 

Jo 
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That is, (f>{L) is the operator with domain 

Dom(0(i)) = S^feL^in) : 

defined by 



l^Xt dEfjiX) < oo 



(2.2) 



(0(i)/,5)L2(a) 



0(A) dE{X)f,g 



0(A) dEf-giX). 



L2(f2) "'0 

These considerations allow us to define the following operators: 

The heat-diffusion semigroup generated by L: with domain X^(J7), 



-tL 



-tx 



dE{X), t > 0. 



We have the contraction property in L^{n): ||e /||^2(q) ^ 
The fractional operators i'^, for < o- < 1: with domain Doni(L'^) C Dom(L), 

The negative povirers i for > 0: 



(2.3) 



L- 



X-" dE{X) 



1 



n^) Jo 



dt 



Proof of Theorem 

1. First we prove that u{-,y) e L'^{Q) and, for aU g e i^(f2). 



(2.4) (u(.,y),.g(.))i. 
For each i? > we let 

UR{x,y) 



1 



r(^) Jo 



(e-*^(L-/),.g)^,(„^e 



dt 



2 dt 



V) Jo 



Since / G Dom(L'^), e~*^(L'^/) e L^(f2). Moreover, e^^r /t^^'^ is integrable near as a function of t. 
Then, using Bochner's Theorem, (|2.2p . the fact that dEf^g{X) is of bounded variation, and the change 
of variables t = r/A, we have 



y2 dt 



r(fT) Jo 

^ r Te-'^A- di?/,(A) e-* 



1 



Via) 
1 

rv) 



oo pR 



^0 



2 ^+ 

e-*\tA)'^e-fe _ di?^^^(A) 



"'0 



dr 



e-W'e-i^^ y dEf^g{X), 



so that 



(?^i?.(-,2/),3(-))L2(n) 



< 



Via) 



oo /-oo 



"'0 



dr 



e-W^ - d\EfJ{X)<\\fh,^^)\\gh.^^y 



Therefore, for each fixed y > 0, uii{-,y) is in i^(i7), and y)llL2(-f2-) < 

The last calculation shows that limi^-^^R^^oo {uR2i-,y) — URii-,y), g{-)) ]^2(^q^ — 0. Then, for any 
sequence of positive numbers, with /• oo, the family {u^^j (•, ?/)}^gpj is a Cauchy sequence 
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of bounded linear operators on L^(il). Thus, uii{-,y) — > u{-,y) weakly in L'^{n), a.s R ^ oo, and 
u(-,y) e L^(ri). Moreover, 



1 r /-^ z dt 



(u(^y),5(0)L2(j^) = ^lini^(ui?X^y),.9(0)L^(o) = ^li^f(^y^ e '^{tXye j dEf^g{\) 

1 r r z d< 1 r r „ ... _z 



1 nOO nOO fif 1 /* C 

1 r,-tL,^-^^^ - 



(e-'^(L'^/),g)i2(a)e 



where the limit can be taken inside the integral because the double integral converges absolutely. 
Hence, (|2.4p follows. 

2. Next we show that u{-,y) E Doni(i), that is, 

lim / £2M:iykli^,<^(.)\ exists for aU g e L^n). 



2 dt 



s-i-0+ \ S 



As e~*^ is self-adjoint, by (|2.4p we have 



r(a) io ^ -"^/L^cn)- ii--- 

Hence, (j2.4L (I2.2p . Fubini's Theorem and dominated convergence give 



L2(f2) r((T) Jo \ S / L2(0) 



dEf.^iX) e-^ 



n^) Jo Jo s ' t'-" 

1 e^'^e'^^X" -e-'^X" dt ,^ 

IV) io io ^ e-«^.i.,,(A) 

1 / / , -tA-" dt 



na)J, J, - ^ ^^^^ ' ti 



3. We check the boundary condition (jl.Sp : for g e ^^(rj), by (12. 4p 



4. The function u is differentiable with respect to y and 



(2.5) uyix,y) = — - / e-'^{L'^f){x) dy{e-^) --^ = —- e-'^{L-f){x) 
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Indeed applying (j2.4p . dominated convergence and Bochner's Theorem we get 

/ u{-,y + h)-u{-) ,A 1 r / -tLtT'y,\ \ ( '^^ 



5. The function u verifies the extension equation (jl.6p . Observe that the integral defining Uy in 
(j2.5D is absolutely convergent as a Bochner integral, and it can be differentiated again with respect 
to y. Hence, 



L2(0) 

^ /'OO 



A- di?/.,(A) 



dt 



II j'\-"-iL' = {Lu{-,y).g{-)) 



6. Let us check (|1.8I) . Note that, for all g e L'^{^), by (|2.4p and the change of variables t = y"^ /{4,r), 

therefore, since limt^o+ *'^^'^/' 5)l2(q) = i-^'^ f ^ 9) (n) ' by dominated convergence, we obtain the 
first identity in (|1.8p . Using (|2.5I) and the same change of variables, the second equality of (|1.8p 
follows analogously. 

7. Finally, we derive the Poisson formula (ll.9p . By (|2.4p . (|2.2p . Fubini's Theorem and the change 
of variables t — j/^/(4rA), we get 



-1 /"CJO /"CJO 2 y^f 



The last equality is due to Bochner's Theorem. 

The second identity of (|1.9p follows from the first one via the change of variables r — y"^ / (4t). □ 

In what follows, we assume that the heat-diffusion semigroup is given by integration against a 
nonnegativc heat kernel Kt{x, z), t > 0, a;, z G f2, that is, 

-tL 



e-*^/(x)= / Kt{x,z)f{z) dr^iz). 
Jn 

Since e~*^ is self-adjoint, Kt{x, z) — Kt{z, x). The second assumption we make is that the heat kernel 
belongs to the domain of L, and dtKt{x, z) = LKt{x, z), the derivative with respect to t is understood 
in the classical sense. This implies that 

dt I Kt{x, z)f{z) dr^{z) = / dtKtix, z)f{z) dTj{z), f e L^n). 
Jn Jn 

Motivated by concrete examples, we add the hypotheses that given x, there exists a constant Cx and 

e > 0, such that \\Kt{x, •)IL2(f^) + \\dtKt{x, •)llL2(n) < C^i^ + t')t~' ■ 
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Theorem 2.1 (Poisson formula). Denote hyVyf{x) the function u{x , y) given in (|1.9p . Then: 
(1) We have Vyf{x) = I Py{x,z)f{z) dr]{z), where the Poisson kernel 

(2.6) p;(^,^):^ Ktix,z)e~'ir 



4'-r(a) Jo ' ' ' 



is, for each fixed 2 e f2, an (Q,) -function that verifies (jl.6p . 

tL , 



(2) sup^>o IT'^-Jl < sup,>o \e~'\f\, m fl. 

(3) If e~*^ has the contraction property in LP(J7), then W'Py fWj^p^^^^) — II'^IIlp(O)' f'^^ '^^^ 2/^0- 

(4) //lim,^o+ e-^'^f = f m LP{n), then lim^^o+ V^f ^ f m LP{n). 

Proof. The integral formula in (1) can be verified by using (jl.9p . Bochner's and Fubini's Theorems. 

In order to see that the Poisson kernel satisfies (|1.6p . we begin by showing that it belongs to the 
domain of L. By the assumptions established on the L -norm of the heat kernel, P^i-,z) G L'^i^l), 
for each z, and, by Bochner's Theorem, 

^"''^^^^•'")=4%)i, ^-''K,{x,z)e~i^ s>0. 

With this, we have 

e-''^P^{;z)~P-{;z) y^- p e-^^Kt{-,z)-Kt{-,z) _Z dr 



We use the Mean Value Theorem, the fact that Kt{-,z) £ Dom(L), and the contraction property of 



e to get 

=\\Le i^t(-,z)||^2(o) = r LKt{-,z)\\^ 

<l|Li^t(.,z)ll^2(^,) = ||a,i^,(.,z)ll^2(o) <c^.(i + ^^r^ 

Hence, the Dominated Convergence Theorem (for Bochner integrals) can be applied in (|2.7p to see 
that the limit as s — )■ 0+ of both sides exists, and 

= a.(if.(x,z))e-^^. 

Now we are in position to check that Py{x,z) verifies p.6p . Note that, by dominated convergence, 
the derivatives with respect to y of Py{x,z) exist and can be computed by differentiation inside the 
integral sign in (j2.6l) . Then, using integration by parts, 

dyPy^ix,z)+dyyPy^X,z) = -^^^ .(x , z) C' ^ " ^) ^ 

2cr roc 2 ^+ 

thus (1) is proved. (2) follows from the second identity of (|1.9p . The contraction property of the heat 
semigroup gives (3): 



I poo ^2 ^ 



Finally, observe that 



\V^f-f\\ < ^ r \\e~^\f - f\\ e-^^ 
\'y-l J \\Lp{n) - Y{a) Jo " ■'\\Lp{nf ^i-<, 



so (4) follows. □ 

Remark 2.2. Note in (jl.9p that, when a = 1/2, V^^^ f = e^*^/ is the classical subordinated 
Poisson semigroup of L acting on / (see [TTl p. 47 and 49]). 
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Proposition 2.3 (Fundamental solution of (|1.6p ). The function 
(2.8) ^:{z,y) = ^J^ Kt{x,z)e''^ 

satisfies equation (|1.6p . = and 

(2-9) lim y'~''^dym;y)J{-)) ^ ^^^m 

Proof. As in the proof of Theorem II .11 it can be checked that for each x, 
' I , . dt 



1 /""^ -tL / N 
r(a) 7o i'-" 
and that (|2.8p satisfies (|1.6I) . Differentiation with respect to y inside the integral in (|2.8p can be 
performed to get 

With this we obtain (1^ : 

-r(i-a) ^ n-a) 



□ 



Remark 2.4. It can also be proved that 



Proposition 2.5. Lef v{x,y) — y^^'^"uy{x,y), where u solves (jl.6p . T/ien v is a solution of the 
following "conjugate equation" 

(2.10) — Lv — - — — Vy + Vyy — 0, in V, X (0, oo). 

Proof. The calculation is analogous to the one given in JJ, with the obvious modifications, and we 
omit it here. □ 

Remark 2.6. As in [2] the fundamental solution (|2.8p and the "conjugate equation" (I2.10p (which 
coincides with the conjugate equation given in [5] when L = —A) can help us to find the Poisson 
kernel (I2.6p . Indeed, we want to write u{x,y) = Vyf{x) — f^Py{x,z)f{z) dri{z) where the Poisson 
kernel Py{x,z) must be a solution of (I1.6P for all z and limj^^o+ T^yfi^) — fi^)- The right choice 
would be 

(2.11) P;(x,z) = r(_^(73|)2(r-a) y''''•''''^9y^i'''iz,y)^C^^,y'-'^'-'^^^y^l-^z,y), 

since it solves the "conjugate equation" (j2.10p with 1 — cr in the place of a (thus it verifies (jl.6p ) and 
by (|2.9p and the choice of Ci-a, 



hm / yi-2(i-)a,vl/i-(z,y)/(z)d,y(z) = /(x). 

A simple calculation shows that (|2.1ip coincides with 
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In the following discussion we shall assume that the operator L can be factorized as L = D*Di, 
where Di — ai{xi)dxi + bi{xi), is a one dimensional (in the ith direction) partial differential operator, 
and D* is the formal adjoint (with respect to drj) of Di . See examples at the end of this Section. In 
this case we give a definition of n conjugate functions related to the Poisson formula for u. 

Let Ecr = —L+ ^y" dy + dyy. Then the factorization 

n 

^ - E + y-^'-^^^dy^y^'^^dyl 

i = l 

suggests the following definition of Cauchy-Riemann equations for a system of functions u, ui, . . . , w. 
17 X (0, cx)) — R such that E^^u — 0: 



(2.12) 



y^-'^"dyU^Dlvi + ---+Dlvn, 

A" = 2/~'-^~^'^''9yDj, i = l,...,n, 

DkVi ^ DiVk, i,k = l,...,n. 



Proposition 2.7. Let u be a solution of E^u = m 57 x (0, oo). //ui, . . . , u„ verify (|2.12p then each 
Vi solves the ith conjugate equation 

(2.13) E{_^Vi = -Lvi + [D*,Di\vi~- — — dyVi + dyyVi = 0, i = l,...,n, 

where [D*, A] = A* A - AA- 
Proof. 

-Lv, + [D*,D,]v, = - ^ DlDkV, - D,D*v, = - ^ DlD.Vk - D,D*v, = -A [Y.^>>A 
= -A (y^-^^dyu) = -y^-^-dy{D,u) = -y^-^^dy{y-^^-^^^dyV,) 

□ 

Remark 2.8. The ith conjugate equation (|2.13p is not the same as the "conjugate equation" (|2.10p . 
They will coincide only when [D*^Di\ ~ 0. This is the case if L = —A: the conjugate equation 
established in [5] is equal to each ith conjugate equation (|2.13p . 

Proposition 2.9. Fix z G f7 and choose u{x,y) = Py{x,z). Then a solution to (|2.12p is given by 
the n conjugate Poisson kernels defined by 

(2.14) v,{x,y):^Ql\x,z)^^-^D,J^ Kt{x,z)e-'- i^l,...,n. 

Proof. From (|2.11l) and the second equation of (I2.12p we have Ci-aDidy'^^jT'^ [z, y) — dyQy''^{x, z), so, 
in view of (j2.8p . Qy'^{x, z) can be chosen as in (I^Ti]) . Clearly AQ;;''(a;,z) = A<3;;^''(a;,z). Moreover, 
the first equation of (j2.12l) also holds: 



DlQl^'{x,z) + --- + DlQl^-{x,z)^^^^Y.^D*D, Kdx,z)e-^ - 



, ^, , / LKt(x,z)e « — 

2 r^r^. , .vl dt 
. ^, , / OtKt(x,z)e «i* — 

' Kt{x, z)e fj ' 



4'^r(cr) Jo ^ ' ' \At J ti+'^ 
:yi-2-9,p;(x,z). 
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□ 

Corollary 2.10. The Poisson integral of /, u{x, y) = Vy f{x), and the n conjugate Poisson integrals 
oj f defined by 

(2.15) v.,{x,y) s Q'^/fix) / Q^;^ {x , z) f {z) dv{z) = ■j:;^D, / e-''^f{x)e-^ - 



4-r(a) 'Jo ' t- 



for i = 1, . . . ,n, solve (|2.12p . 



Remark 2.11. When a = 1/2, Q^^^'^ f{x) is the ith conjugate function of / associated to L, see [TT| 
and [12]. A natural question arises: what is the limit of Qy''''f{x) as ?/ — 0^7 The answer is contained 
in the next result. 

Theorem 2.12. For each x £ Q, 

Proof. From the expression of Qy'''f{x) in (j2.15D and (I2.3L 

□ 

Remark 2.13. The conclusion of Theorem 12.121 can also be obtained from the following observation: 
except for a multiplicative constant, the last formula of (j2.15l) is just the D^-derivative of the solution 
of the extension problem (jl.6p for L^~'^ with boundary value L^'^^^'^^ f{x) (see (jl.7p ). For a — 1/2, 
Theorem 12. 121 establishes the boundary convergence to the Riesz transforms DiL^^/^ (which in case 
L = —A are the classical Riesz transforms S^;. (— A)^^/^). See |Tl] and [I2] . 

Examples 2.14. We present now some examples of operators L for which our results apply. 
The Laplacian in R": Observe that, when f G S, 

(2.16) e'^fix) = f Wt{x - z)f{z) dz, Wt^x) = — i— e"^. 

The Poisson formula given in 2 is recovered: use the change of variables = r, in 

to see that the Poisson kernel in this case is 



p.,-A, y'" re-^^ir^ dt _rin/2 + a) y^" 

r [X, z) 



4'^r(cr) Jo (47ri)"/2 ti+'^ 7r"/2r((T) / . 



(^\x - zf' + 2/2^ 



The function Py^'^' ^{x,z) is the classical Poisson kernel for the harmonic extension of a 
function to the upper half space. The Cauchy-Riemann equations read 



' y^-^'^dyu = - (a,,fi + • • • + d^^vn) , 



(2.17) 



dxiU ^ y^'^^^^'^^dyVi, i = l,...,n, 

dxk'Vi = dxiVk, i,k = l,...,n. 

The case a = 1/2 is the classical Cauchy-Riemann system for the n conjugate harmonic 
functions to u. In dimension one (|2.17p reduces to 



dxu - y-'^^-^^'^dyi 



which already appeared in [7]. 
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Classical expansions: L can be each one of the operators arising in orthogonal expansions, 
like the Ornstein-Uhlenbeck operator (Hermite polynomials and Gaussian measure dri{x) = 
e-l"l' dx), 

-A + 2X-V = J2 + 2a;^) ; 

i 

the harmonic oscillator (Hermite functions and Lebesgue measure dr^{x) = dx), 
-A + Ixl' ^\Y1 [("^-' + + ^'') + + + ^^)] ' 

i 

the Laguerre operator (Laguerre polynomials and measure dri{x) — Y[i^'i^^~^^ ) 

Jacobi and ultraspherical on (—1, 1); etc. 

We would like to point out that in these cases, due to the existence of smooth eigenf unctions, 
the proof of Theorem \l.l\ can be performed as an exercise of convergence of orthogonal systems, 
and it makes it technically simpler. 
Elliptic operators: Let L be a positive self-adjoint linear elliptic partial differential operator 
on L^(fl), with Dirichlet boundary conditions, and bounded measurable coefficients. Then the 
heat kernel exists, and it verifies our assumptions stated before Theorem \2.1i Even more, its 
heat kernel has Gaussian bounds ^ p. 89] . We can also consider Schrodinger operators with 
nonnegative potentials in a large class [21 Section 4.5]. 

3. Existence and uniqueness results for the extension problem 

In this section we derive tlie concrete solution of ttie extension problem in the case of discrete 
spectrum. We also find solutions with null Neumann condition. This is done by using classical 
Fourier's method. 

Let {(pkJkeng be an orthonormal basis of L^(rj) such that L<j)k = Xk'Pk, Afe > 0. Recall the definition 
of L°' given in (fTTO)) . 

3.1. theory. Let / e L'^{fl) and look for solutions u to (|1.5p - (|1.6p of the form 

(3.1) u{x,y) = ^Ck{y)(l)k{x). 

k 

Then for each fc > we have to solve the following ordinary differential equation: 

-XkCk + - — — 4 + c'l^O, y> 0, 

y 

with initial condition Cfc(O) — (/, (pk)- According to [51 p. 106], this last equation has a general solution 
of the form 

(3.2) Cfc(y) -y-Z,(±zAy'2/), 

where Zg- is a linear combination of Besscl functions of order a. To have uniqueness of the solu- 
tion include the boundary condition limy^oo u{x,y) ~ 0, weakly in L'^lTt), which translates to the 
coefficients as 

(3.3) lim ck{y) = 0. 
From [51 p. 104], can be written as 

(3.4) Z,{z) = AiJ,(z) + A2iJ«(z) = B,Mz) + B2Hi^\z) 

(3.5) = CiMz) + C2J-AZ) = DiH^^\z) + D2H^''\z), 
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where J^- denotes the Bessel function of the first kind and Ha^^ and -ff^'^' are the Hankel functions. 
To fulfill condition p.3p we need to review the asymptotic behavior of the Bessel functions. When 

|arg z| < TT — (5, 

1/2 



(3.6) 



TTZ J 
TTZ 



2cr7r + TT 

4 

2(T7r + TT 



l + 0{\z\-')) 



\ 8z 



+ oi\zr) 



'-^ {l + 0{\z\-')), 
-) [l + 0{\z\-')) 



Note that for purely imaginary z, Jaiz) — ^ oo exponentially and h'^\z) — >■ or c» depending on the 
sign of the imaginary part of z. Putting z = i'^k'^V in p.2l) we see that the only possible choice as 



solution is the first linear combination of p.4p as soon as Ai = 0: Ck{y) = A2,ky'^ H^\i\^J'^y). If K^r 
denotes the modified Bessel function of the third kind then Ha^\iz) = i-K^^i^"^^ Kc{z) and 



To determine A2^k use the initial condition. The asymptotic behavior of K(j{z) as z — > reads 
(3.7) 



i^.(z)«r(a)2--i^. 



So that, when y 0, Ck{y) w ^2,fe2'"7r-ir'"-ir((T)A, 



-(t/2 



Therefore 



2,fc 



,-l+cr 



2'^r(cr) 



Thus 
(3.8) 

Since as Izl 



1/2 



{l + Oi\z\-') 



the series in p.ip . with as in p.Sp . converges in L^(r2) for each y e (0, oo). Finally, p.7p implies 
that (jl.Sp is fulfilled in the i^(ri) sense. 

On the other hand, by using the properties of the derivatives of (see 6j p. 110]) and (13.71) . as 
2/ — >■ we have 



dix'J'y) 



iXl/'yrKAX'J'y) 



d{\Ty) 



1 1 9^21-'^ 

—y 

2a 



{f,^k)i-m'J'yrK.-i{\'J'y)\'J' 



2-" 



-aT{a) 
2-2'^r(l - a 

-'^r(a) 



(/,0,)Afe/'Af yi-ifi_.(Af y) 
{f,^k)Xl^^\'jV-''ni-<7)2 



(Afy)- 
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As a consequence, 

the limit taken in the L^(ri)-sense (see (jl.lOp ). 

3.2. Local Neumann solutions. Let us find a solution to (jl.6p such that 

(3.9) -!- lim y^-^''uy{x,y) = 0, for all x e n. 

2(7 J/— >o+ 

Writing u{x,y) ~ J2k dk{y)4>k{x), condition p.9p implies that limj^_5.o+ y^^^'^ d'f.{y) ~ 0. Therefore as 
(see (IX^ ) 

d'.iy) = (zAfc)!-'^— ^ \i^Xl/'yrZ,i^\l/'y)] = zAf y-Z._i(zAf y), 



we require 

y'-^^d'.iy) = iXl^V-'^Z^^^itXl^'y) 0, y ^ 0. 

When z ^ (see [B]), 

2'^i (1 + Cr) iTT Z'^ ITT Z'^ 



Then, as y — ^ 0, 



TT 



TT 



but 

/•\1/2m- 

l_o- ^ /•\l/2 \ ^ l*^fe ^ 2-2cr , n 



Consequently, choose the first linear combination in (13. 5p with Ci = 0. Thus 

d/c(y) = C2,ky'^J-aii>^l^^y), 

verifies liuiy^Q y^~^'^ d'i.{y) — 0. So u formally reads 

y)^y"^ C2,kJ-a{iXl/'^y)<j)k{x). 

k 

In order to have a convergent series (at least for small y), let us determine C2,k- Taking into account 

dnSl) it is enough to fix i? > and put C2,k = Ce"^^^'^. 

In this way we obtained a solution u to equation (jl.6p in O x (0, i?) that satisfies the required 
property p.Op . 

4. The Harnack's inequality for ff'" 



To prove the Harnack's inequality in Theorem 11.21 we first study the problem (jl.5p - (|1.6p for the 
harmonic oscillator L = H = —A + |a;| posed in M" with the Lebesgue measure drj = dx. 
Fix 1 < p < oo and > 0. Define the space 

(4.1) L^;, = <( u:M"^R: ||u|l^P = f / ^^^^^ dz\ < oo 



Then L% CS'. 
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The heat semigroup generated by H (see [12] ) can be given as an integral operator 

/• /• — [^^ |a: — 2:|^ coth 2i+a:-2: tanh 

(4.2) e-*«/(x) = / Gt{x,z)f{z) dz = / . f{z) dz. 



(27rsinh2i)"/2 

We collect some useful facts about e^*^ in the next Proposition, whose proof is postponed to the end 
of Section [S] 

-tH - 

up c 



Proposition 4.1. For f E L^, the heat semigroup e *^/(x) is well defined and 



(4.3) \e-'"f{x)\ < C'- '^'J, X G R", t > 0, 

where p > depends on p and N. Moreover, [dt + H)e^*^ f{x) = for all x G M" and t > and /or 
= 1, . . . 

(4.4) |9..(e--^/)(x)|<C ^ ' ' |a...,(e-^/)(x)|<C-' 



^(n+l)/2 ' rai.XjA J A J\-^ ^(n+2)/2 

// / is aZso a function in some open subset O C K" then limt_j.o e~*^/(a;) = f{x) for all x E O. 

In the particular case we are considering in this Section, Theorem 11.11 takes the following form, in 
which the relevant observation is that all identities are classical. 

Theorem 4.2. If f E L^j^ is a function in some open subset O C M" then 

2a poo 2 

(4.5) u{x,y):^J— e-'"f{x)e-- 



4'^r((T) Jo ' ' 

is well defined for all x E K", y > 0, and 

—HxU + - — — Uy + Uyy = 0, 171 M" X (O, oo ) ; 

lim u{x,y) = f{x), forxEO. 

In addition, for all x E O, 

1 1 r°° df 

(4.6) — lim y^-^"uy{x,y) = —— / {e~'" f {x) ^ f {x)) 

^ ^ 2CTa^o+ ' AT (a) J a ^ ^ ^ t^+'' 

Proof. Estimate (j4.3l) implies that the integral defining u is absolutely convergent and Uy and Uyy can 
be computed by taking the derivatives inside the integral sign. Moreover, by using (|4.4I) . we have 

,,2cr /-oo 2 rl-f 



HM^.y) = -TTW^, He-'" fix) 



in the classical sense. Hence, for each x E M", u verifies the extension problem in the classical sense. 
To check that (j4.6p is classical, we begin by recalling that 



9 / .,2 



Thus 



2t) 



2a" 2aA''r{a) Jo \ 2t J t^+<' 

As we shall see later (Remark l5.12l) the integral in (|4.6p is absolutely convergent for all x E O and / 
as in the hypotheses. Therefore (|4.6p follows. □ 

Remark 4.3. In Section[5]we will see that for f E L^j^ H C^{0), H"^ f is weh defined and 

1 , ^ dt 
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Remark 4.4. Theorem 14.21 is valid if H is replaced by —A and the function /, with the same 
smoothness in O, belongs to L^r := ^'^'^ discussion on (— A)"" given in Section[5l 

Lemma 4.5 (Reflection extension). Fix R> and xa E K". Let u be a solution of 

-HxU + ^ Uy + Uyy = 0, 171 M" X (O, R) , 

with 

(4.7) lim y^~'^" Uy[x,y) = 0, for every x such that \x — xa\ < R. 

Then the extension to R" x (— i?, R) defined by 

u{x,y), y>0; 
u{x,-y), y < 0; 



(4.8) u{x,y)-^ 
verifies the degenerate Schrddinger equation 

(4.9) divi\y\'-'^Vu)-\y\'-"'\x\'u = 0, 



in the weak sense in B 



\{x,y) e M"+i : \x~xo\^ < R^}. 



Proof. A nontrivial solution u can be found with the method of Subsection l3.2l since the eigenfunctions 
of the harmonic oscillator H are the Hermite functions ha, a E Nq , with corresponding eigenvalues 
Xa = 2\a\ + n (see Section[S|). Given ip e C^{B) we want to prove that 



/ := / fvu • Vip + |a;|^ uip] \y\^ "^^ dx dy = 0. 

J B 



For 5 > we have 

/ = 



divdyl""^ "^"^ cpWu) dx dy + / f Vu • V(/3 + |a;|^ j lyl'^ ^'^ dx dy 

Bn{\y\>s} JbdUvKS} ^ ' 

— / i-pb^^'^" Uy(x,S) dx ^ / [y u ■ \I 1^ + \x\^ ui~p\ \y\^~'^" dx dy. 

JBn{\y\=s} JbdUvKS} ^ ' 

As (5 — ?► 0, the first term above goes to zero because of (|4.7I) and the second term goes to zero because 
^|V'u|'^ + \xf' lyP"^'^ is a locally integrable function. □ 



Proof of Theorem M .B. Let u be as in Theorem 14.21 Since / is a nonnegative function, from (j4.2p and 
(|4.5I) we see that u > 0. Because of Remark 14.31 its reflection (|4.8p satisfles Lemma 14.51 Note that 
(|4.9I) is a degenerate Schrodinger equation with A2 weight w = and potential V = |?/|^~^'^ \x'^ 

such that V jw e £^ locally for p large enough. So we can apply the result of 0] to obtain the 
Harnack's inequality for u and thus for /. □ 

5. POINTWISE FORMULA FOR AND SOME OF ITS CONSEQUENCES 

The semigroup language adopted in Section [21 allows us to get the exact pointwise formula for the 
fractional Laplacian (— A)'^ on M". The constants involved in the definition are computed exactly in 
an easy way. 

Lemma 5.1. For f £ S, 

Proof. The first identity follows by Fourier transform. From the fact that e*^l(a;) = 1 we can write 

r°° df r°° r dt 

(5.2) / (e'^fix) - fix)) —= IWtix- z){f{z) - f{x)) dz — = h + h, 
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where 

h r / Wdx - y){J{z) - f{x)) dz 

I — 1^ 

and Wt is the heat kernel for the Laplacian (j2.16p . Use the change of variables s = ^ to see that 

1 dt _ 4^(11/2 + cr) 1 



^^■^^ /o (47ri)«/2 " " ti+- W2 |x-zr+'"' 

So, since / is bounded, Ii converges absolutely. Passing to polar coordinates, 

By Taylor's Theorem, l^,,^{f{x + rz') - f{x)) dS{z') = Cnr^Af{x) + 0{r^), thus 



1 /-oo „-^- J, /■I 



I/2I < C„,A/(:r) ^^75:^ T '^'^ " CnAf(x).<y dr = C„,a/(x),^, 

and I2 converges. Therefore apply Fubini's Theorem in (|5.2p and (|5.3p to get (|5.ip . □ 

Remark 5.2. Lemma [5T] gives the exact value of the positive constant Cn.a in (|1.2p . Observe that 

, , -4'"r(n/2 + CT) , 

(5.4) c„,<, = -75-7 ^0, as cr 0+ or cr ^ 1 . 

7r"/^l (— cr) 

When / e 5 it is clear (by Fourier transform) that linio-^i- (— A)'^/ —A/. The next Proposition 
shows that this is in fact valid for / G C^. Note that if / € 5 then, from (|TTT|) . {-A^f ^ S, but still 
{-AYf e C°°. It can be checked that for every /3 € NJf the function (1 + |a;|"+^'")i:>'^(-A)'^/(x) is 
bounded. Therefore the set ■= |m : R" — > M : — /u„ i^zY'+'i" < oo| (which is 

in (|4.ip '). consists of all locally integrable tempered distributions u for which {~A)°'u can be defined. 
If / e Lct is in an open set O then it can be proved that {—A)'^f is a continuous function in O 
and its values are given by the second integral in (jS.ip . For all the details see |9j and (lOj . 

Proposition 5.3. Let f £ C'^{B2{x)) D L°°(M") for some x eW\ Then 

lim {-Arf{x) = ~Af{x). 

cr— i-l- 

Proof. Fix an arbitrary e > 0. Since / G C'^{B2{x)) there exists (5 = (5^ > such that 



(5.5) l-D /(w) - /(w')| < for all w,w' e Bi(x) such that \w-w'\<5. 

Write (-A)-/(a;) = c„,,(/ + //) where / - /|^_^|^, g!^^."] d^. We have |/| < C„a-i<5-2- 

so that from (|5.4p . c„ ,j/ as cr ^ 1^. Using polar coordinates, Taylor's Theorem and recalling 

//= / r-^-^" [ {f{x)- f{x-rz')) dS{z') dr 



J\z'\ = l 

r-^-^" [ Rif{x,rz')dS{z')dr 



s 



-A/(.T)(n/2 + l)7r"/2r2 /" / 



2r(n/2 + 2) 7,,,, 1 V ' ' ' 2 



i?i/(a;,rz') " ^(^V(a^)^',^') H5(z') 



dr 



' / f i?i/(a;, rz') - '^{D^f[x)z\ z') ] dS{z')dr 



-A/(x)(n/2 + l)7r"/252-2'^ /- / / „ ,^ r 



4r(n/2 + 2)(l-a) ^ y,,,,^! V ^ ' ' 2 
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where Rif{x,rz') is the Taylor's remainder of first order. Then (|5.4p entails 

-Af{x)ain/2 + l)r{n/2 + a)S^-^- ^ ^ , K2 + l)r(n/2 + 1) 
^"-^^^ = 4i-r(n/2 + 2)r(2-a) ^ -^-^(^^ r(,V2 + 2) = 



as cr — > 1 . Finally, by (jS.Sp . 



Rifix,rz')- i^{D^f{x)z',z') < C^r^e and I//2I < C„<52-2'^(l-CT)^ie. 
Therefore Umo._i.i-i |c„,cr//2| < C„£. □ 

Remark 5.4. For / e C^(i32(a;)) H Lo- the second identity in (|5.ip is valid (the idea is to use the 
continuity of f as in the proof of Proposition [221). 

We shall now discuss the definition of the fractional harmonic oscillator H"' and the pointwise 
formula for H°'f{x). The eigenfunctions of H (see [12]) are the multi-dimensional Hermite functions 
defined on K" as ha{x) = ^aix)e~^^^ a € Nq, where are the multi-dimensional Hermite 
polynomials, and Hha = (2 |a| + n)ha. Note that ha G S. The set of Hermite functions forms an 
orthonormal basis of L^(R"). Let f E S. The Hermite series expansion of / given by 

CO 

(5.6) ^(/, ha)K ^Y.Y.^^^ '^")^- 

a k=0\a\=k 

with {f,ha) = J^,t fha dx (which converges to / in L^), converges uniformly in R" to /. This 
uniform convergence is a consequence of the fact that ||^a||2^oo(Rn) < C for all a € Ng and the 
following estimate: for every m £ N, 



{2\a\ -l-n)™ ~ {2\a\ +71)™' 
since H is & symmetric operator. If / G 5 then 

(5.8) e-'"f{x) = Y,e-'^^\^\+-\f,K)ha{x), t > 0, 

the series converging uniformly in E". By the given estimates on H/JqU^oo and |(/, the series 
defining the fractional Hermite operator 

(5.9) i7"/ = ^(2|a|-t-n)'^(/,/i„)/i„ 

converges uniformly in M". 
Lemma 5.5. For f £ S, 

Proof. Let Cq = (/, ha). Because of the uniform convergence of the series of (|5.6p . (|5.8p and (I5.9P we 
get 



(e--/(x) - /(.)) ^ 



[ (^^e-(2|"l+")c„/^.(x)-^c„/^.(.)^ ^ 



dt 
t^ 

= T{-a)Y,{2\a\+nYcaha{x) = r(-a)i?V(a:)- 

a 

□ 

We have the following important Lemma whose technical proof is given at the end of this section. 
Lemma 5.6. H'^ is a continuous operator on S . 
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Lemma 15.61 together with the symmetry of H"^ on S (that can be easily verified via Hermite series 
expansions) allow us to give a distributional definition of Bf^: for u E S' , define H'^u E S' through 

Therefore H'^ is well defined for all functions u that are tempered distributions. In particular, u can 
be taken from the space of (j4.ip . 1 <p<oo, iV>0. 

Recall the expression of Gt given in (14. 2p and the fact that (see [5,) 

(5.10) e-'"l{x) = - L__ e-^l-l' < 1. 

^ ' ^ ' (cosh2i)"/2 

Define the nonnegative functions 

1 f°° dt 1 f°° dt 

(5.11) FAx,z):=^^J^ G,{x,z)^, B^(^) := __ (e-*^l(^) _ l) _. 

Theorem 5.7. Let f be a function in that is C^{0) for some open subset O C M". Then f 
is a continuous function in O and 

H''f{x)^SJ{x) + f[x)B,{x), xeO, 

where 

(5.12) S„f{x)^ f FAx,z){f{x)- fiz)) dz. 

In (|5.12p we see that H°' is a nonlocal operator. Before giving the proof of Theorem l5.7l we establish 
some easy consequences. 

Theorem 5.8 (Maximum principle for H"^). Let f be a function in L^ that is in an open set 
O C K". Assume that / > and f{xo) = for some xq G O. Then H'^f{xo) < 0. Moreover, 
H'^ f{xo) — only when / = 0. 

Proof. By Theorem EUl since f,F^>0, 

H^fixo) = / (/(.To) - f{z))F,{xo, z) dz + f{x^)B„{xo) = f f{z)F,{xo, z) dz < 0. 

If f{z) > in some set of positive measure, then the last inequality is strict. □ 

Corollary 5.9 (Comparison principle for H'^). Let f,g E i^v ^ C^iO) be such that f > g and 
f{xo) ~ g{xa) at some xq E O. Then H'^f{xQ) < H'^g{xo). Moreover, f(xt)) = H'^g{xo) only 
when f = g. 

We devote the rest of this paper to the proofs of Lemma [5761 Theorem [5771 Proposition [43] and to 
complete the missing details at the end of Section [31 

Proof of Lemma 1 5. 61 Define the first order partial differential operators 

. _ 9 A - ^ -1 

Ai : — — h Xi, A—i : — -;- \- Xi, i — 1, . . . , n. 

oxi oxi 

It is well known that 

(5.13) A,ha,{x) = (2a,)'/'/ia-e.(a;), A_,/i„(a;) = (2a, + 2f/^ho.+,^{x), 

where is the ith coordinate vector in Nq (see [I2|)- This implies that H'^ f E C°° and for all fc G N, 

(5.14) A,,---A,^H''f{x)^Y.^2\a\+nY{f,h^)A,,---A.Mx), n ^ ±\, I = I, . . . ,k, 

the series converging uniformly on M". Since 

Ai + A—i Ai A—i d . 

2 2 "a^' « = 
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for each multi-index 7, /3 € Nn we can write x'^D^ = xl^ ■ ■ ■ xZ" — „F^'^^ r as a finite linear combina- 

tion of operators Ai and A^i. Therefore, to check that x^D^H"' f £ L°° , it is enough to verify that 
for each fc e N, A,^ ■ ■ ■ Ai^H'^f e L°° where {ii, . . . , ik} C {-1, 1}. The identities in (|5.13p easily 
imply the following commutation relations for Hcrmitc functions and thus for f G S: 

AiH'^J = {H + 2YAif., l<i<n\ 
A.H''! = {H -2YA,J, ~n<i<-l. 

Here {H ± 2)M,/ := ^^(2 |a| + n ± 2)-(AJ, K)K. Hence, in ([531, 

A,,--- A,,H^f - ^(2 |a| + n + 2jr{g, h^)h^, 



for some j £1^ and g := Ai-^^ ■ ■ ■ Ai^f G S. For m G N sufficiently large, as in (j5.7p . we have 
^(2|a| +n + 2jr(g,/i„)/i„(x) 



Therefore x^D^^Wf G L°°. Moreover, 

< C (seminorms in S of [Ai^ ■ ■ ■ Ai,, f))~C (seminorms in 5 of /) . 



□ 



For the proof of Theorem 15.71 we need some estimates on Gt , F„ and . First we derive some 
equivalent formulas for these kernels. Consider the change of parameters due to S. Meda 

(5.15) t = = ifogi±£, tG (0,00), sG (0,1), 

z 1 — s 

that produces 

(5.16) -^ = d/v(s):= -j^, t £ {0,00), s£ [0,1). 

(l-.2)(ilogi±|) 

Then the heat kernel in (|4.2p can be written as 

/I 2 \ "/2 

G.(.)(x,z)=(i^] e-M^I^+^l^+il-^n, s£iOA), 
and, from (|5.11l) and (|5.16p . 

Lemma 5.10. For all s £ (0, 1) and x,z £ R", 

(5.17) G,(,)(a;,z)<C^ e-^^e-^. 



/n particular, 

(5.18) Gt(,)(a;, z) < (1 - sr/\--^e--^e- — 

\x — z\ 



Proof. The second estimate in the statement follows immediately from (j5.17p . Note that 

G,(,){x,z) < C ( i^y^^'e-^e-M^I^+^l^+ll-^n < C ( V^' g-^e-^l^-^H^+^l. 



We prove the second inequality above. Assume first that |x — z| < |a; + z\. Then by minimizing the 
function e[s) := f |a; + z\^ + |^ k - z\^ for s £ (0, 1) we get g-sH^+^l'+^l^-^ll < ^-^,\x-z\\x+z\ _ 
the case \x + z\ < \x - z\ we have e"4H^+^l'+7k-^l1 < g-iVl^-^l' = f.-h\x-z\\x-z\ < ^-^l\x^z\\x+z\ 



20 



P. R. STINGA AND J. L. TORREA 



for all s e (0, 1). To obtain estimate (|5.17p proceed as follows: if a; • z > then |a; + z| > which 
gives 8 < e^s'^H^^^I; if a; • z < then |a; — z| > \x\ and in this situation 



e s« e 



e 16s e 16s < e 16= e 



□ 



Observe in (15.16^ that 



(5.19) 



ds 



, S 0, dfla{s) 



ds 



(l-.s)(-log(l-s)) 



l + CT 



s - 1. 



Lemma 5.11. For all x,z E 
(5.20) F^{x,z)< 



C 



1=^11 



|a; — z\ 

Moreover, B„ e C°°(E"). 
Proof. Estimate (|5.18p gives 



:n+2cr 



e ^' e ^' and B, 



.{x) <c(l + \xf 



Fa{x,z) < C 



e ^ 



\x - z\ 



(l-s)"/2e-^ d^,(s). 



Then (|5.19D implies 

»l/2 



(1 - s)"/2e-"'c^ d/i^(s) < C / e"^ 



1/2 



if \x- z\< 1; 



and 



1/2 



/fi la; — z| la; — z| 

(1 - s)"/2e-^^ d/v(s) < 



ds 



1/2 (1-S)(-I0g(l-S))l + 



e if |a; — z| > 1; 
Ce 



thus the first inequality in (|5.20p follows. 
Apply in ((glI7| to obtain 



(5.21) 



^^"l{x) 



1 - s' 
1 + s2 



i/2 



e 1+^ 



Then, up to the factor _p^^_^^ , we can write 



B,[x) 



l-s' 

1 + S2 



2 \ "/2 



e 1+= 



d^i,{s) + / (e^lf:^!"!' - l) d^i,(s) ^I + II. 



To estimate / and // we use (|5.19l) and the Mean Value Theorem. That is, 



\I\<C [ 
Jo 



1/2 



1-S^ 

1 + s2 



2\ »/2 



^ + /' df,„{s) <C [ 

Jl/2 Jo 



,1+cr 



c = c. 



For // we consider two cases. Assume first that |a;| < 2. Then 



|//| < C 



1/2 



e 1+= 



ds 

b1+ct 



1/2 



1/2 



dHais)<C I \x\^s^ + C<C. 



In the case |x| > 2, 



C|a;|'" + C 



ds 

bI+o- 



+ / rfM^(s) < 



1 (l-s)(-log(l-s)) 



1+cr 



- fog 1 - 



<c\ 
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since — log(l — s) ^ s as s — > 0. Therefore the second fact of (|5.20l) foUows. is differentiable since 
the gradient of the integrand in its definition is bounded by 

„2\ "/2 



21x1 



1 + V 1 + s 



1-s' 



e-TTI^I^I <C\x\seL\{0,l)■,d^l,is)), 



thus we can differentiate inside the integrah 



Vi..(.)=2.y^ ^(^) e-TT.N^.M.). 
For higher order derivatives we can proceed similarly. 

Proof of Theorem \5.1\ Take first f G S. Since e~*^f (x) is not a constant function we write 



□ 



^ - /(.)) ^ 



Gt{x,z)f{z)dz-f{x)] — 



dt 



Gtix, z){f{z) - fix)) dz + fix) / Gtix, z)dz-l 



dt 



df r°° df 

Gtix, z)ifiz) - fix)) dz— + fix) I (e-*^f (x) - 1) - 



Gt{s)ix, z)ifiz) - fix)) dz dfi^is) + fix)B„ix). 

Due to Lemma [?751 the first integral above is well defined and converges absolutely. Write the integral 
in the last line as Is + Is<^ with Igc = J\x-z\>s Gt{s)ix, z)ifiz) — fix)) dz d^a-is), for some S > 
(in this step d is arbitrary, but we will fix it later). Estimate (15.201) implies that /^c is absolutely 
convergent and {Is^l < G ||/|lioo(j{n'). Pass to polar coordinates in Ig: 

,2 \ "/2 



Is 



Airs 



.ih\x+z\'' + i:\x-z\^ 



ifiz) - fix)) dz dfiais) 



\x-z\<S 



Vf-s2 



2 \ n/2 „5 



Airs 



\z'\ = l 



'M-+"'\\fix + rz') - fix)) dSiz')drdfiAs). 



To estimate /gn-i :— J\^n^i e * 1"^^+^^ \ ifix + rz') — fix)) dSiz') use the Taylor expansions of / and 



f\z'\ 

H™l and cancel out terms: 



\z'\-. 



(e-tl2-l' +i?oV^(a;,rz')) (V/(x)(rz') + r/)) dSiz') 



\z'\ = l 



e-tl2-l'i?i/(x, rz') + R^^six, rz')V/(x)(rz') + R^^^ix, rz')Rifix, rz') dSiz'). 



Since \Ro^psix,rz')\ < s^/^r and \Rifix,rz')\ < ||^V|Loo(b^(^)) we have |/s"-i| < Cr^. Thus 



\Is\< 



1 _ "/2 



Atts 



< c 



1 r 1 



„n/2 



r"-ie-5^ |/5„-i| drdfi^is) 
s 

r.n+1 



e dfi„is) dr < G r 



1 



dr = GS^ 



r 



Hence Is converges. The conclusion follows, for / G 5, by Fubini's Theorem. 

Now assume that / G L^, 1 < p < oo, N > 0, is a G^ function in O. Then H"' f is well defined 
as a tempered distribution. Fix an arbitrary x G O and take (5 > so that Bsix) C O. Observe 
that the integral in (|5.12l) is well defined: just apply Taylor's Theorem (as above) in Is and the 
condition together with (|5.20p in Igc. Let e > 0. There exists a sequence fk G C^(M") such that 
ll^^/fe|L~(B,(x)) - \\^'^A\l=-(Bs{x)) ^' h converges uniformly to / in Bsix) and fk ^ / in 
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the norm of L^, as fc — > cxi (use moUifiers and multiplication by a smooth cutoff function). Since B„ 
is a continuous function, JkB^r converges uniformly to fB„ on Bs{x). Let < p < S/2 he such that 
for all k 



F<y{x, z){Jk{x) ~ fk{z)) dz 



< |, and 



F,{x,z)ifix)- fiz)) dz 



Bp(x) 



e 

<3- 



For k sufficiently large, by Holder's inequality. 



F,{x, z){h{x) - h{z)) dz - / z)(/(x) - f{z)) dy 

B-(x) JB-ix) 



< 



<\fk{x)- f{x)\ F,ix,z)dz+ F^{x,z)\.fk{z)~ f{z)\ dz 

Jb-{x) Jb-{x) 



< 



C{\fk{x)-f{x)\ + \\h-fh.^) 



< 



Thus 



Safk{x) 



F„{x,z){f{x)- f{z)) dz 



in Bs{x). But H'^ fk — > H"^ f in S' . By uniqueness of the limits, S„f{x) coincides with the integral in 
(j5.12l) . Moreover, H'^ f is continuous in Bs{x) because it is the uniform limit of continuous functions. 

□ 



Proof of Proposition pTT] By (|5.15p . (|5.17p and Holder's inequality, 

C 11/11.., 



-t{s 



< 



,n/2 



{1 + Izi'fP' dz 



^ C — 



n/2 



For note that if < s < i, then s < t{s) < |s. The equality dte~*"f{x) = dtGt{x, z)f{z) dz 
is valid if the last integral is absolutely convergent for all t in some interval. But dtGt{x, z)f{z) = 
—HxGt{x,z)f{z), therefore we have to verify that /jj„ HxGt(s){x, z)f{z) dz converges absolutely for 
all s in some interval. This last statement is true since 



|Va:Gt(s)(x,z)| 



< 



5I/2 



-c[s\x+z\'' + i\x-z\^] 



and 



\DlGt(s){x, z)\ 



< 



1-s' 



1 \ n/2 



\ ^-c\s\x+z\^ + \\x-z\^ 



which give estimates similar to (I5.17P for VGt(s) and D'^Gt(^s)- Hence dtc = —H^e 

and (|4.4I) follows. Observe that t{s) — > if and only if s 0. For a; e O we have 



e-*(^)^/(x) - fix) 



< 



Gt(s){x,z){f{z) - f{x)) dz 



\f{x)\ e-*(^)^l(x) - 1 



The last term above tends to as t{s) — > because of (|5.2ip . Let 5 > be such that Bs{x) C O. 
Then, as / e C\Bs{x)), 

Gtis){x,z){f{z)- f{x)) dz <C 



lBs(x) 



e cs 



\z — x\ dz < C 



Bs{x) 



e Cs 



Bs{x) \Z - X\ 



— c?z — 0, 
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when s — > 0, by the Dominated Convergence Theorem. On the other hand, 



Gt{s){x,z){f{z) - /(x)) dz 



< 



i/p 



e 



c 

=: / X //. 

Clearly I < oo and, by dominated convergence. 



{i + \zn^p {i + \z\')Np 



dz 



^e-'^il + lzlY"' dz 



i/p' 



1/p' 



IKC 



-e -^{l + lzffP dz] ^0, ass^O. 



= {x) \x-z\ 



□ 



Remark 5.12. If / e n (7^(0) then, for each x eO, 



e-'" fix) -fix) 



dt 



e-*'^-^'>"f{x)- fix) d^iAs)< 



Indeed, by gj 



-tH 



dt 



fix)-fix)\ -^<Cix) 



1 dt 



5 log 3 



< OO, 



and 



5 1°S 3 fjf 



1/2 



< C 



1/2 



Gt(s)(2;, z)(/(z) - fix))dz 



ds 



IT^ + C|/(x)| 



1/2 



1 - 

1 + s"^ 



2\ "/2 



e 1+ 



ds 



Both integrals above are finite: the first one by the arguments in the proof of Theorem 15.71 (Taylor's 
Theorem) and the second one because of the Mean Value Theorem. 
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